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Abstract 

We study whether the pseudoscalar glueball mass in full QCD can differ from the 
prediction of quenched lattice calculations. Using properties of the correlator of the 
vacuum topological susceptibility we derive an expression for the upper bound on the 
QCD glueball mass. We show that the QCD pseudoscalar glueball is lighter than the 
pure Yang-Mills theory glueball studied in quenched lattice calculations. The mass 
difference between those two states is of order of l/N^.. The value calculated for the 
^ QCD glueball mass can not be reconciled with any physical state observed so far 
in the corresponding channel. The glueball decay constant and its production rate 
in J/ip radiative decays are calculated. The production rate is large enough to be 
studied experimentally. 



PACS numbers: 12.60. Jv; 14.80.Ly; 11.55.Hx. 
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Introduction 



Glueballs are one of the intriguing theoretical predictions of QCD [|T|. The search 
for these composites is a long-standing problem of theory and experiment. We now 
have detailed experimental studies of the resonances in the mass region up to 2.3 
GeV, as well as great progress in lattice QCD calculations. A number of interesting 
particles have been detected 0. Some of them, the ?7(1410, 0""*"), the /o(1500, 0"*"+), 
the /o(1710, 0^"*") and the ('(2230,2"*""'') appear to have a rich gluon content. (For a 
recent discussion of the phenomenology of these composites and a full set of references 
see [^). The glueball candidates can be compared to lattice QCD predictions[|. These 
calculations argue in favor of the following hierarchy of glueball masses: the 0+''' 
glueball is the lightest one with mass about 1.5 — 1.7 GeV [|], ]10[; the 2+"^ state. 



having mass 2 — 2.2 GeV P], is the next one; finally, the 0^''' pseudoscalar glueball, 
being predicted in the lattice calculations to have mass 2.3±0.2 GeV is the heaviest 
one. Various studies of pure Yang-Mills (YM) theory also support the picture outlined 



above. In ref. |]TT| the theorem was proved that in pure gluodynamics the pseudoscalar 
glueball is heavier than the scalar one. Instanton calculations [|1^] also confirm this 
picture. 

Notice that all the theoretical facts listed above are firmly established results of 
pure YM theory. One might wonder whether this picture is affected when quarks are 
also included in the theory. Recent analyses flOl, llll, [§], |l^ of the scalar glueball 



candidates indicate an important mixing with the nearby qq resonances. This leads 
to a modification of the mass spectrum and the decay constants of the glueball states 
0, Ijl^. The mass shifts due to mixing are approximately 100 MeV or so, and the 



lattice predictions are in good agreement at this level with the experimental data for 
the 0++ and the 2++ channels [§, . 



However, the situation for the pseudoscalar channel is problematic. One can ex- 
pect that, because of the axial anomaly, quarks are crucial for ^ channel physics 



ISf As we mentioned above, the lightest "•" glueball predicted by the lattice calcu- 



lations has mass about 2.3 GeV [Hi. On the other hand, there is evidence that the 



The earliest theoretical predictions were based on QCD sum rule calculations (^-||^. 



77(1410) can have a substantial gluonic component 0] Q Hence, an important question 
is how the quark degrees of freedom may shift the glueball mass in the ^ channel 
and whether one can identify the QCD glueball with any observed state. We examine 
here these questions. 

In low energy QCD the ?7'(958) makes the dominant contribution to correlators 
in the ^ flavor singlet channel. It is well known that the mass and decays of the Vj 
meson are strongly affected by the gluonic sector of the theory. The axial anomaly 
and instantons play a key role in generating the mass and decay constant of the Vj 
(see also the papers ^T|, ||22[ and refs. therein). This suggests that the pure 



+ glueball should be also affected by the singlet qq admixture in full QCD . 



In what follows we will refer to this mixed glueball-gg state as a QCD glueball in 
distinction with the pure YM glueball and the physical v( . 

We are going to derive an inequality between the mass of the QCD pseudoscalar 
glueball and the mass of the YM glueball as measured in quenched lattice calculations 
0. In the derivation we closely follow the argumentation of Witten p3[ and Veneziano 
but keep truck of the finite iVc effects by using the method of QCD sum rules 0]. 
The decay constant and mass for the QCD pseudoscalar glueball will be determined. 
We show that the QCD glueball is lighter then the one of pure YM theory. Finally, 
we will look at J j'^ radiative decays and predict the production rate for the QCD 
glueball state. 

1. QCD Sum Rules and the Witten- Veneziano Relation 

In this section we are going to study the properties of the correlator of the vacuum 
topological susceptibility. In Minkowski space-time this is defined as 

X(g' = -Q') = ( J^) ^ / e^^^(0|TG;,G;,(x)G^^G^^(0)|0)rf^x, (1) 

where = l^^i/a/jGjJ,^ and i§2^G^yG^y is the renormalized composite operator with 

g being the running coupling constant. 

^Note that structure originally identified as 77(1440) is in fact two states, the 77(1410) and the 
7;(1490), [|l^-[^. The lighter r;(1410) seems to have a rich gluon content 1^, while the heavier 
77(1490) is dominantly ss state possibly with some admixture of glue 



This correlator has a different behavior in the vicinity of = depending on 
whether it is evaluated in full QCD or in pure Yang-Mills theory p3|. If light quarks 



are included in the theory, Xqcd(O) is proportional to the product of the quark masses 
and vanishes in the chiral limit. This fact is related to the absence of the theta angle 
in massless QCD. In general, the theta angle can always be rotated away by an 
appropriate chiral transformation of a massless fermionic field. On the other hand, 
there are no massless fermions in pure YM theory. As a result, the explicit theta 
dependence cannot be removed ||2^ and Xym(O) turns out to be a nonzero number. 
Let us consider the dispersion relation for the function x{Q'^)/Q'^ 

X(0) - X{Q^) r P{s)ds 

Q2 Trio s{s + Q^) 

Following the standard QCD sum rule approach Q], the spectral density for this 

correlator, p, can be decomposed into two parts. The first one consists of resonance 

(pole) contributions and the second one is determined by the perturbative expansion 



/■°° p{s)ds , . 

/ ~/ 7^ + subtractions. (2) 

Jo sis + Q-^) 



where the CnS are resonance residues, m^'s are corresponding masses, and sq denotes 
the continuum threshold; p is given by the perturbative expansion of the correspond- 
ing correlator. We are going to work in the next-to-leading order (NLO) of the per- 
turbative expansion. In the case at hand, the leading contribution to p is scale and 
scheme independent. However, the next-to- leading term depends on the renormaliza- 
tion scheme. In leading order p is fixed by the diagram which contains only gluon 
propagators. This diagram is the same for QCD and pure Yang-Mills theory. Hence, 
in lowest order Pqcd = Pym- However, this relation does not hold in the NLO. There 
are quark loop corrections to the gluon propagator in QCD. Hence, the result for p 
in NLO QCD differs from the one of pure YM theory by quark loop contributions. 
In order to fix the scale/scheme ambiguity of the NLO corrections it is convenient to 
adopt the Brodsky-Lepage-Mackenzie (BLM) scale fixing procedure . In the BLM 



scheme NLO quark loop insertions into the gluon propagator are summed up into the 
redefinition of the effective scale of the strong coupling constant. Hence, in the BLM 
scheme, the perturbative expansion for Pqcd in the next-to-leading order formally 



coincides with that for pym- Keeping in mind that we have adopted the BLM scheme 
one can write down the following relation for the spectral densities in QCD and pure 
Yang-Mills theory ^ 

S-a + 5^)s'^as'. (3) 

Before we turn to the resonance part of the spectral density let us make an important 
comment. It deals with the definition of the running coupling constant in QCD and in 
YM theory. The expression for as{s/fi^) depends on the number of flavors Nf present 
in the theory. In pure YM theory Nf = and the coupling constant of this theory 
differs from the one defined in full QCD. However, our goal is to stay maximally close 
to what is used in quenched lattice calculations. In those calculations only gluon 
degrees of freedom are taken into account. But this is not the whole story. Quenched 
lattice calculations effectively include some of the virtual quark effects through the 
formal substitution of the QCD running coupling with Nf = 3 instead of the coupling 
of YM theory with Nf = (see ref. for this discussion). We are using this formal 
method through the paper. In particular, the theory to which we refer as pure YM is 
actually the theory with some of the virtual quark loops effectively included through 
the use of the QCD running coupling constant instead of the pure YM running 
coupling. Thus, YM theory in our context refers to the model which has the full 
QCD coupling constant, but nevertheless, differs from true QCD by the absence of 
qq bound states and the absence of quark condensate effects. Let us stress again that 
these conventions differ form the normally used ones (adopted for example in ref. 
p4| ) and are motivated by the quenched lattice calculation procedure. 

After this remark let us turn to the resonance part of the spectral density. This 
part for YM theory is assumed to be saturated by the pure glueball state Gq, and for 
QCD by the r]' meson and QCD glueball state G. The expressions for the spectral 
densities are 

Pym(s) = fco^Go^is - m^J + pym(s)^(s - Sq), 

■^In this case the BLM scheme leads to the better convergence of the perturbation expansion. 
The NLO corrections are large in the AIS scheme 0. 



Pqcd — Pym 
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Pqcd(s) = flm%5{s - my) + f^,m^,S{s - mj,) + Pqcd{s)9{s - si), (4) 

where Sq and Si denote the continuum thresholds for YM theory and full QCD re- 
spectively. Other quantities in these equations are defined as follows 

{0\^,G%G%W) = (5) 

with rriG being the QCD glueball mass, m^o being the YM glueball mass, and fc and 
fco the corresponding decay constant^ 

Before we turn to the application of the QCD sum rule method let us first compare 
the Operator Product Expansions (OPE) for the quantity B^Q"^) = — 
QCD and in YM theory. In leading order only the gluon fields contribute in both 
cases. The results of calculation of these OPE's can be found in refs. 0, 0, @: 



where 



D'^ = D^ + Oia,m,{qq)), D'^ = D, + 0{al{qq)^). (7) 

The instanton contributions in eq. (6) are suppressed as where n ~ 12 

Since for practical calculations we use eq. (3), the NLO perturbative corrections are 

not explicitly shown in eq. (6) for brevity. Let us make some comments about the 

quantity -Bqcd- mentioned above, the perturbative part of this correlator in the 

NLO is the same as that of By^. However, there are nonperturbative contributions 

in -Bqcd('5^) which do not appear in the expression for i?YM^(<5^)- Those are related 

^The gluonic operator -^f^p^-GG appearing in eq. (5) has an anomalous dimension so that the 
constants /gq , fc ^^^id are not renormalization group invariant quantities p6| . For any of three 
/'s one can construct the renormahzation group invariant constant by means of a finite multiphcative 
renormahzation of /'s E^, p8[. 
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to the quark condensate. The first such contribution modifies the term in 

the OPE for -Bqcd- The new contribution is proportional to agrriq^qq). The next 
nonperturbative correction related to the quark condensate appears at order and 
yields a new term in addition to the operator Oq. The additional term is proportional 

to a^'^ 



Now we can turn to the QCD sum rule analysis. It is useful to apply the Borel 
transformation (with the Borel parameter denoted by M^) to eq. (2) with the 
phenomenological part on its l.h.s. and the OPE on its r.h.s. This leads to the 
following sum rules 



oo 



s Jo s ySzTT'^y V 

+Xym(0), 

io s Jo s \ ^ M2 J 

Taking now the limit — > oo we get[] 

- Xym(O) + foo^k = ^ + - ( J^) ' D„ (8) 

/^-^^ + /X' = ^ + -(j^)'^4- (9) 

We have mentioned already that our goal is to study how the QCD glueball mass 

differs from the one of quenched lattice gauge theory. However, eqs. (8,9) alone 

are not yet enough to determine whether that difference really exists. We need an 

additional relation. One way to get the new relation is to use the dispersion relation 

for the correlator of the topological susceptibility xiQ"^) itself 

1 r°° p{s)ds 
IT Jo s + Q"^ 

In order to get rid of the subtractions let us use again the Borel transformation. 
Applying this transformation to the correlator one gets the following relation 



, 1 1"°° pis)ds 
Y(Q I = — / — ■ + subtractions. 

IT Jo S + Q2 



TT JO 



^In this way one derives the finite energy sum rules ||29|] modified by condensate contributions 



where we dropped the subscripts distinguishing QCD from YM theory. Taking the 
hmit ^ oo we obtain the following relations 

/l.< = ^-''(3f^)'c«. /G"i + /K' = ^-''(3f^)'«;.- (10) 

Having eqs. (8-10), one can use them to find the relations between the quantities 
defined in YM theory and in full QCD. For that goal we are going to use the Witten- 
Veneziano arguments, but along with leading terms we keep corrections of order 1/A^c- 
In the large Nc limit fc'rrvQ —>■ /gq'^Go' ^^"^ accordance with the Witten-Veneziano 
(WV) relation [^], p4|p] %ym(0) = —f^im^,. As a result the left hand sides of eqs. 



(8) and (9) are equal in the large Nc limit. On the other hand, the second terms on 
the r.h.s. of eqs. (8,9) are also equal in that limit (difference between them is of order 
1/Nl). So, we come to the conclusion that lim7Vc^oo(si — sq) = 0. The same result 
can be drawn from consideration of eqs. (10) in the large Nc limit. 

Since we are going to keep 1/Nc corrections as well as leading order terms let us 
introduce the following parametrization for the continuum thresholds 

8 , , 

si = so + — , (11) 

where b is an unknown quantity which is of order of the unity in the large Nc limit. 
In what follows we are going to solve the system of equations (8-11) keeping track of 
'VjNc corrections. It is convenient to introduce the following notations 

^4(6) = vr J-^) ^4(6), ^4(6) ^^(j;^) ^4(6)- 

Based on the rules of the large Nc expansion one can see that 
D',-D,^Oim„—), D',-D,^OC^). 



"in our conventions x is defined in Minkowski space-time and turns out to be a negative quan- 
tity. One can turn to Euclidean space-time in eq. (1) making the substitutions xq —ix^ and 
GG — !■ iGG. The Euchdean quantity, which is used in the lattice calculations, is a positive number 
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Using the notations given above, and neglecting all contributions of order l/N"^ and 
higher, in the chiral limit one can rewrite the system of eqs. (8-10) as follows 

-Xym(0) + /^Xo = ^ + ^4, /^o< = ^-^6, (12) 

/^-^^ + /.X' = ^+^4, /S-^a + /X' = ^ - ^6. (13) 

Before we go further let us make some comments. Two equations given in (12) contain 
only two unknowns /gq and sq (assuming that the YM glueball mass is known 
from lattice calculations 0). Hence, those two equations can in general be solved 
and the values of /g„ and Sq can uniquely be determined. This is done in the next 
section. On the other hand, the two equations in (13) contain three unknowns, fc 
and Si and itlg- So one can not determine uniquely. The only thing one can do is 
to calculate the decay constant fc and mass ma for chosen values of the continuum 
threshold Si which are dictated by previous analyses of the flavor singlet pseudoscalar 
channel. This calculation is also carried out below. Before we turn to the numerical 
simulations we can try to extract some analytic relations for the QCD glueball mass 
and decay constant studying the system of equations (12-13). The relation between 
the continuum threshold parameters (11) allows one to set the equations (12-13) as 
following 

flml + f^,ml, = -Xym(O) + /^„m^„ + ^.o, (14) 

Now the key observation is that in pure YM theory there are not light meson states, 
thus the decay Gq Stt does not occur in this theory. On the other hand, this decay 
should easily go in QCD. Hence, the continuum threshold for QCD si should be less 
or equal to the continuum threshold for YM theory^ 

si < So or 5< 0. (16) 

Thus, up to order and in the chiral limit the following inequalities exist between 
the quantities defined in pure YM theory and in full QCD 

fl^l + fn'^ l' < -Xym(O) + fa,ml^, (17) 

^ I am grateful to Glennys Farrar for bringing this line of arguments to my attention. 
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fomh + f^,m^, < fc,m%^. (18) 

Let us point out that eq. (17) yields the Witten-Veneziano relation in the limit of 
infinite Nc- Indeed, in the large Nc limit fc'^G ~^ /go^^Go' ■^o -^i- Hence, in that 
limit eq. (17) is saturated and it turns into the Witten-Veneziano (WV) relation |2£ 



H Xym(O) = -f^^m^,. 



Let us rewrite eq. (18) in the following form 



This inequality allows one to calculate the upper bound on the pseudoscalar glueball 
mass in QCD. The inequality shows that if the values for fc^ and fc are sufficiently 
close to each other, then the QCD glueball is lighter than the glueball of YM theory. 
In the next section, based on numerical studies, we demonstrate that this indeed is 
the case. 

It is interesting to check the large A''^ behavior of eq. (19). Notice that, since 
the anomaly disappears in the limit when ^ oo, there is no flavor singlet meson 
-glueball mixing term anymore in the effective Lagrangian in that limit. One 



should therefore expect to have equal masses for the glueballs in QCD and pure YM 
theory when A^^ oo. Recalling that /^/ oc y/Nc, m^, oc l/Nc, fcj^'ho ~ 1 ^-^d 
substituting these into eq. (19) we get 

Go G 

which is consistent with ones expectation. 

3. Some Estimates and Predictions 

Let us now turn to numerical estimates. First of all let us list all the approxima- 
tions we made deriving eqs. (12-15). There are scheme dependent NLO perturbative 
corrections involved in the derivation. Besides that, we worked in the chiral limit 
neglecting m, d and, most importantly, s quark masses. Hence, the rj — rj' mixing and 
all other nonsinglet pseudoscalar mesons are also neglected [|. Below we show that 
*A11 these corrections are expected to be of order 0{mu.d,s/ ^giuebaii), where Agiuebaii — 1-5 ~ 
2 GeV is an effective scale above which the existence of glueballs become important for hadron 
physics. 
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the contributions of the 77(1295), the 77(1410) and the 7/(1490) in the sum rules are 
rather small and can also be neglected. 

If one knew all the numerical values for the quantities on the r.h.s of eq. (19), one 
would be able to predict the upper bound on the QCD glueball mass. Unfortunately, 
fco and fc are not known. We can use however the sum rules derived in the previous 
section to calculate the value for fc^. Indeed, consider the equations given in (12). 
One can solve this system with fco and sq treated as unknowns. The numerical 
solution for that system yields the following result = (27 ± 3) MeV and Sq = 
7.4 ± 0.5 GeV^. 

In calculating the numbers given above we have used the following numerical 
values for = (29±3) MeV,f}m^, = (957.77±0.14) MeV, = (2.3±0.2) GeV §, 
as(2 GeV) = 0.33±0.05. There are a number of estimates for the gluon condensate in 
the literature (see refs. [Q, |38|). We take the world average value of these calculations 
{^Gl^) = (2.5 ± 0.9) 10-2 GeV^ The corresponding values for and Dq are 
D4 = (4.0 ± 1.7)10-^ GeV^ = (0.7 ± 0.3)10"^ GeV^ 

There are also several lattice estimates for the topological susceptibility Xym(O) 
(for a recent review see |^). For our estimates we take the result of refs. ||33[| , 
Xym'(O) = (175±5 MeV)^, which is also in good agreement with an earlier theoretical 
estimate |4|Q 

Now we can make the prediction for the matrix elements of the gluonic operator 
acting on the pseudoscalar pure YM glueball state with mass uigo = (2.3 ± 0.2) GeV 

(0|^'G;:,G^JGo) = (45±9)GeV3. (20) 

The uncertainty in this result dominantly comes from the error bars associated 

with the value for the topological susceptibility and also with the value of frj' 0. 

tWhich corresponds to Fq = (105 ± 12) MeV take n in conventional normalization for the singlet 
axial current. 

**In order to study whether our results are sensitive to the numerical value of xym(O) we varied 

the value of this quantity from (100 MeV)^ to (190 MeV)'*. Outside of the interval given by these 

numbers the system of equations does not have a solution. On the other hand, inside the interval 

the results are rather insensitive to the value of xym(O), varying xym(O), sq changes in the region 

6.8 - 8 GeV and fco varies in the interval 24 - 34 MeV. 

"''The only available lattice calculations for the matrix elements of the gluonic operators were 
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The same matrix element for the Vj meson state has the following numerical value 
(Ol^^G^^^G^JV) = (8.4 ± 0.8) GeV^ The values of the matrix element for the QCD 
glueball state G will also be given below. Notice that all the numerical results for 
the scale dependent quantities such as the decay constants /gq, /g and the matrix 
element in eq. (20) are to be taken at the normalization point approximately equal 
to the value of the glueball mass. Also, all the results presented above should be 
given 10-20 percent systematic error bars associated with the approximations and 
the method we have used. 

As we mentioned already we checked whether the presence of the ?7(1295), the 
?7(1410) and the ?7(1490) mesons could affect our results. We have included the con- 
tributions of these resonances in the sum rules used above. In order to determine the 
decay constants of these resonances we used the experimental data for the production 
rate of these states in J j^^ radiative decays []TB|, [^. The values of the decay con- 
stants (defined as in eq. (5)) are rather small numbers: fnixm^) = (7.5 ± 1.9) MeV, 
/j7(i4io) = (4.8 — 6.7) MeV and /,,(i295) = (4.5 ± 1.0) MeV. Including these numbers 
into the sum rules one can see that the final results for the glueball mass and decay 
constants change unsubstantially. 

Now let us turn to the calculation of the mass and decay constant of the QCD 
pseudoscalar glueball. Unfortunately, our method does not allow one to determine 
these values uniquely. The set of two equations given in (13) contains three unknowns, 
/g, rnc and s\. Below we present the results of numerical solutions of these equations 
for different values of the continuum threshold s\. The values for s\ are chosen 
between the upper bound determined as Si < Sq = 7.4 GeV^ and the lower bound 
given by Si > 3 GeV^ (below this value the continuum threshold comes very close to 
the resonance mass square and applicability of the sum rule approach breaks down). 
The results are summarized in Table iQ. 



presented in ref. pSl]. However, the predictions of ref. |3q] are in conflict with the theoretical 



estimates in the case of 0++ channel (see the discussion in ref. |12 ) and also with our results for 
^ channel. 

^^^^In Table 1 and also Table 2 below only the error bars associated with the method of numerical 
calculations are given. 
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Si GeV^ 


fn MeV 




si GeV^ 


fn MeV 


TD.r^ GpV 

/ / ''(-(- ^— * ^ V 


7.4 


29 ± 2 


2.27 ± 0.04 


5.0 


21.5 ± 2.5 


1.9 ± 0.05 


7.0 


28 ±2 


2.2 ±0.04 


4.5 


20 ±3 


1.8 ±0.1 


6.5 


26.5 ± 1.5 


2.15 ±0.05 


4.0 


18 ±3.5 


1.73 ±0.12 


6.0 


25 ±2 


2.07 ±0.05 


3.5 


16 ±4 


1.61 ±0.14 


5.5 


23 ±2 


1.97 ±0.05 


3.0 


13.5 ±5.5 


1.47 ±0.20 



Table 1: Sum rule results 



As we see from the table the QCD glueball is lighter than the glueball of YM 
theory. However, it is hard to identify the QCD glueball with the 77(1410). The very 
low value for the continuum threshold is needed in order to have QCD glueball mass 
at about 1.4 GeV. 

Though we can not determine uniquely the value for the continuum threshold from 
our consideration, one can use the estimate given in ref. @] si ~ lOm^ ~ 6 GeV^. If 



this number is accepted, then in accordance with Table 1 the QCD glueball mass is 
niG = (2.07±0.05±0.3 (syst.)) GeV and decay constant fc = (25±2±4 (syst.)) MeV. 

Having these numbers at hand one can predict the J/ip decay width into a glue- 
ball state and photon. J/ip radiative decays are very effective tools in studying the 
spectroscopy on light mesons. In the present case we are going to deal with the pro- 
cesses like J/ip -—>■ R{0^^) 7, where R stands for the resonance being considered. The 
theory of these decays was worked out in ref. One can consider the ratio 

T{J/i; ^ G 7) 



r = 



T{J/ij^7]' 7)- 

This ratio is independent of the J/ip meson wave function and is completely defined 
by the properties of the pseudoscalar mesons produced in the decay. Assuming that 
the decay dominantly goes through the exchange of the intermediate gluons in the 



pseudoscalar state, the ratio r can be rewritten as follows 



§See also the approach developed in H). 
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fn MeV 




k GeV^ 


T ratio 


r keV 


29 ± 2 


2 27 ± 04 


47.2 ± 5.0 


4.2 ± 1.1 


1.57 ± 0.4 


28 ± 2 


2 2 ± 04 


42.8 ± 4.7 


4.2 ± 1.0 


1 57 ± 37 


25 ±2 


2.07 ±0.05 


33.8 ±4.5 


3.7 ± 1.0 


1.38 ±0.37 


21.5 ±2.5 


1.9 ±0.06 


24.5 ±4.6 


2.8 ± 1.1 


1.0 ±0.4 


18 ± 3.5 


1.73 ±0.12 


17.0 ±5.2 


1.8 ± 1.2 


0.67 ±0.45 


13.5 ±5.5 


1.47 ±0.20 


9.2 ±7.5 


0.75 ± 1.4 


0.28 ±0.53 



Table 2: Some predictions 



Using the numerical results listed above in the table one can calculate the ratio r, 
matrix element k = {Olg^G^^^G^^JG) and decay width r(J/V^ ^ G 7) for the QCD 
glueball state. These results are summarized in Table 2. 

The numerical values for the decay width of the J/^ meson into the QCD glue- 
ball and photon are substantially large and this decay can be observed in recent 
experimental studies. 

Conclusions 

Let us summarize our results. We studied whether the pseudoscalar glueball mass 
in full QCD differs from the one determined in the quenched lattice calculations. 

An inequality which sets the upper bound on the mass of the pseudoscalar glueball 
in QCD (eq. (19)) is derived. In order to calculate that bound numerically one needs 
to know the decay constant of the QCD 0"+ glueball and also the mass and decay 
constant of the pure YM glueball. 

The decay constants in this work are calculated using the QCD sum rule approach, 
while the value for the YM glueball mass is taken from lattice calculations. The value 
calculated for the decay constant of pure YM glueball is shown to be important for 
selfconsistency checks of lattice results. 

We found numerically that the mass of the QCD glueball is less than the mass of 
the glueball of pure YM theory. The values for the ^ QCD glueball mass and decay 
constant (for the phenomenologically preferred value of the continuum threshold pa- 
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rameter) are: mc = (2.07±0.05±0.3 (syst.)) GeV and fa = (25±2±4 {syst.)) MeV. 
If these numbers are accepted, then there is no particle discovered so far which might 
be identified with the QCD pseudoscalar gluebalL Further experimental searches in 
the 2 GeV region are needed. In this case the status of the ri{1410) is unclear. 

In order to help resolve this question, we predict the production rate of the QCD 
glueball in the radiative decay of the J/ip meson, r(J/\E' G'-f). For a 2.07 GeV 
glueball r(J/\E' — > G'-f) is about three or four times greater than for the rj' meson. 
Thus, the prediction for the branching ratio for that process is large enough to be 
studied experimentally. 

I am grateful to Glennys R. Farrar for bringing this problem to my attention and 
for valuable discussions and suggestions. The author wishes to thank H. Neuberger 
and T. DeGrand for useful discussions. 

Note added 



After this work was done we become aware of the papers |^| where the QCD 



sum rule method was used to calculate glueball masses and decay constants in scalar, 
pseudoscalar and tensor channels. The calculations in refs. are done (without 
referring to lattice results) using the optimization procedure with respect to Borel 
parameters and continuum thresholds. Our results for the mass and decay constant 
of the pseudoscalar glueball are in good agreement with predictions of EOl. We are 



grateful to M. Schwetz for bringing refs. |40] to our attention. 
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